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We present results of three SLD analyses: our final determination of the rate of gluon splitting 
into bb, an improved measurement of the inclusive b quark fragmentation function in Z° 
decays, and a preliminary first measurement of the energy correlation between the two leading 
B hadrons in Z° decays. Our results are obtained using hadronic Z° decays produced in 
e + e~ annihilations at the Stanford Linear Collider (SLC) between 1996 and 1998 and collected 
in the SLC Large Detector (SLD). In this period, we used an upgraded vertex detector with 
wide acceptance and excellent impact parameter resolution, thus improving considerably our 
tagging capability for low-energy B hadrons. 

1 Introduction 

The vertex representing a gluon splitting into a bb pair, g^bb, is a fundamental component 
of QCD, but the contribution of this vertex to physical processes is poorly known, both theo- 
retically and experimentally. In high-energy e + e~ annihilation the leading-order process con- 
taining this vertex is e + e~ ^qqg^qqbb; information on g^bb can thus be obtained by studying 
e + e~ -^hadrons events comprising four jets, with two of the jets identified as b or b. We de- 
fine the rate g bb as the fraction of e + e~~ — diadrons events in which a gluon splits into bb. The 
measurement of g bb is difficult experimentally since the rate is intrinsically low and the back- 
grounds from Z° —>bb events are two orders of magnitude larger. In addition, the B hadrons 
from g-^bb have relatively low energy and short flight distance and are difficult to identify using 
standard tagging techniques. We present a world's best measurement of g bb using our vertexing 
ability to find low energy B hadrons. 

According to the factorization theorem, the heavy quark fragmentation function can be 
described as a convolution of perturbative and non-perturbative effects. For the b quark, the 
perturbative calculation is in principle understood, and the non-perturbative effects have been 
parametrized in both model-dependent, and model-independent approachesB The most sensitive 
experimental determination of the fragmentation function is expected to come from a precise 
determination of the B hadron energy distribution. We have previously reported the first mea- 
surement of this distribution over the entire allowed kinematic range; i here we update that 
result and exclude many models and place stringent bounds on the form of the distribution. 

Extending the fragmentation analysis to look at the energy correlation between the two B 
hadrons in a Z° —>bb event father constrains theoretical work. It can provide a direct test of the 
QCD factorisation theorem^ and further probe phenomological models. For example, models 
that accurately predict the inclusive energy distribution might not find the correct correlation, 
and vice versa. We present an analysis that measures the correlation and compares it with a 
leading order QCD calculation. 



2 Selecting B Decay Vertices 



A general description of the SLD detector can befound elsewhere, as can the trigger information 
and selection criteria for Z° -^hadron eventsoa The B sample is selected using topologicallv 
reconstructed secondary vertices based on the detection and measurement of charged tracksQ 
To reconstruct the secondary vertices, the space points where track density functions overlap 
were found in three dimensions. Only vertices that are significantly displaced from the IP were 
considered to be a possible i?-hadron. With SLD's good track impact parameter resolution 
and precisely know interaction point, the vertex finding gives excellent 6-tagging efficiency and 
purity. In particular, the efficiency is good even at low 5-hadron energies, which is especially 
important for detecting g— and for shape resolution of the fragmentation function. The mass 
of the reconstructed vertex, M c ^, is calculated by assigning each track associated with it, the 
charged-pion mass. Because of the tiny SLC IP error and the excellent vertex resolution, the 
transverse momentum Pt of tracks associated with the vertex relative to the line joining the 
IP and the secondary vertex is well-measured. The mass of the missing particles can then be 
partially compensated by using P t to form the "Pt corrected mass", Mp t = J M^ h + P 2 + \P t \. 



3 The Rate of Gluon Splitting into bb Pairs 

Candidate events containing a gluon splitting into a bb pair, — > qqg — > qqbb, where the 
initial qq can be any flavor, are required to have 4 jets (Durham algorithm, y cu t = 0.005). A 
secondary vertex is required in each of the two jets with the smallest opening angle in the event, 
yielding 1514 events. This sample is dominated by background (S/N ~ 1/10), primarily from 
Z° — > bbg(g) events, and events with a gluon splitting into a cc pair. 

In order to improve the signal/background ratio, a neural network technique is used. Nine 
observables are chosen for inputs to the neural network. The inputs include the the Pp-corrected 
masses of the vertices and the angle between the vertex axes: b jets have higher P^-corrected 
mass than c/uds jets; many Z° — > bb background events have one 6-jet which was artificially 
split by the jet-finder into 2 jets that tend to be collinear. The neural network is trained using 
Monte Carlo samples of 1800k Z — > qq events, 1200k Z — > bb events, 780k Z — > cc events and 50k 
g — > bb events. A cut on the neural network output of > 0.7 keeps 79 events in the data, with an 
estimated background of 41.9 events. Using this and the estimated efficiency for selecting g — > bb 
splittings of 4.99% yields a measured fraction of hadronic events containing such a splitting of: 

g h l = (2.44 ± 0.59(stat.) ± 0M(syst.)) x 10~ 3 (1) 

The result is consistent with previous results and the QCD prediction! of 1.75 x 10~ 3 . 



4 The Inclusive b Quark Fragmentation Function 

Candidate events for this analysis were found by requiring a secondary vertex in an event 
hemisphere with each hemisphere of the event being treated independently. Placing the condition 
that the vertex must have an M p t > 2 GeV and a distance from the LP. > 0.1 yields a B hadron 
sample of 42,093 from the data with an estimated purity of 98.2%. Assuming all the tracks 
associated with the vertex are pions we calculate the charged track energy, E c h, and the total 
momentum vector P of the B vertex; we must then find the energy of particles missing from the 
vertex. We asssume PTmissing = Prcharged and Mb = 5.26GeV leading to an upper bound on 
the missing mass M : M$ max = M% - 2M B ^M 2 ch + P t 2 + M 2 ch . The true missing mass M$ ue 
is often rather close to Mo max , and we use the latter as an estimate of M§ ue . We can then 
solve for the longitudinal momentum of the missing particles using kinematic considerations, 



and hence the missing energy from the vertex, Eq. The B hadron energy is then Eb = E c h + Eo. 
Since < M^™ 6 < Mo max , the B energy is well-constrained when MQ max is small. For this 
reason, in order to obtain good energy resolution, we cut events from our sample using a value 
of Mg m(M , that depends on Eb to achieve a B selection efficiency nearly independant of xb- A 
total of 4,164 vertices were found after all cut criteria, with an estimated purity of 99.0%. The 
full kinematic range is covered by these events. 

After background subtraction, the distribution of the reconstructed scaled B hadron energy 
is compared with several heavy quark fragmentation models. Within the context of the JETSET 
Monte Carlo the Bowler, Lund, and Kartvelishvili models are consistent with the data as is the 
stand-alone UCLA generator. We also compared the data with a set of ad hoc functional forms 
of the xb distribution in order to estimate the variation in the shape of the xb distribution. 
Of those we tested, only the Peterson function, two ad hoc generalizations of the Peterson 
function and an 8th-order polynomial are consistent with the datall We use the eight consistent 
fragmentation models and functional forms to unfold the raw data distribution for selection 
efficiency and bin migrations. Their average is taken as the central value in each bin and the 
RMS value as the model dependent error. The fully corrected data are shown in Figure [l|. 

3.5 i 1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 1: The corrected single inclusive weakly deacying B hadron energy distribtion. The points represent the 

data, and the envelope is our total error. 

5 The Two Dimensional BB Energy Distribution and Correlation 

The previous analysis has placed considerable constraints on the form of the b fragmentation 
function, and on phenomological models. SLD, however, can go further and look at the cor- 
relation between the 2 leading B hadrons. Events are selected by requiring that exactly two 
secondary vertices are found in different jets (Durham algorithm, y cu t = 0.015). Either vertex 
is required to have an M p t > 2 GeV and the same (or otherwise) vertex, is required to be 
seperated from the LP. by > 0.1 cm. An upper limit on Mq of the remaining sample ensures 
an energy resolution better than 20%. A cut on the angular separation between the two decay 
vertices avoids the problem that sometimes a genuine vertex is split into two separate candi- 
dates. After all selection criteria we have 19283 events in the data in which both the B and B 
are reconstructed with an estimated purity and efficiency of 99.6% and 36.6% respectively. We 



group these events in bins of cos<j) where is the angle between the two vertex axes. By taking 
moments of the 2-dimensional energy distribution we have a measure of the correlation between 
the two hadrons. Within each bin we calculate the double moments Dij(i,j = 1,2,3). We use 
our Monte Carlo to correct the data for selection efficiencies and detector effects. In order to 
compare with the QCD calculation we take a ratio and calculate normalise the Dij to D\\. The 
ratios of these moments, are compared with a leading order QCD calculation in Figure 0. 
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Figure 2: Ratio of the double moments Pij/Pn of the 2-dimensional B hadron energy distribution. Data are 
shown as points with stat. error bars only. A leading order QCD calculation of the ratios are shown as solid lines. 



6 Conclusion 



We have presented the results of three analyses, all of which make use of SLD's exceptional 
B-tagging performance. Our final determination of the rate of gluon splitting into bb pairs per 
hadronic Z° decay, g bb , is a world's best: g bb = (2.44 ± 0.59(stat.) ± 0.34(syst.)) x 1CT 3 . We 
have updated our measurement of the inclusive weakly decaying B hadron energy distribution 
and exclude many models of b quark fragmentation. We obtain a precise measurement of the 
average scaled B hadron energy: < xb >= 0.710 ± 0.003(stat) ±0.005(syst) ±0.004(modeZ). We 
have a first measurement of the energy correlation between the leading B hadrons in Z° decays, 
the results of which are consistent with a leading order QCD calculation. 
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